function r= psnr1(Io,Im)
imagen=Io;
marcada=Im;



if (size(imagen)~=size(marcada)) 
   error('Las imagenes deben tener el mismo tamaño') 
end 
if (~isrgb(imagen) & ~isrgb(marcada)) 
   [m n] = size(imagen); 
   A=double(imagen); 
   B=double(marcada); 
   sumaDif=0; 
   maxI=m*n*max(max(A.^2)); 
   
   for u=1:m 
      for v=1:n 
         sumaDif = sumaDif + (B(u,v)-A(u,v))^2; 
      end 
   end 
   
   if (sumaDif==0)  
      sumaDif=1; 
   end 
   S=(maxI/sumaDif); 
   S=10*log10(S);	 
else  
   error('Las imagenes no pueden ser en color') 
end 
r=S;
end